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Abstract. We study the ground state of the model Hamiltonian of the trimerized S = 
1/2 quantum Heisenberg chain Cu3Cl6(H20)2 • 2II8C4SO2 in which the non-magnetic 
ground state is observed recently. This model consists of stacked trimers and has three 
kinds of coupling constants between spins; the intra-trimer coupling constant Ji and 
the inter-trimer coupling constants J2 and J3. All of these constants are assumed to 
be antiferromagnetic. By use of the analytical method and physical considerations, 
we show that there are three phases on the J2 — J3 plane (J2 = J2/J1, J3 = J3/J1), 
the dimer phase, the spin fluid phase and the ferrimagnetic phase. The dimer phase 
is caused by the frustration effect. In the dimer phase, there exists the excitation 
gap between the two-fold degenerate ground state and the first excited state, which 
explains the non-magnetic ground state observed in Cu3Cl6(H20)2 ■ 2H8C4SO2. We 
also obtain the phase diagram on the J2 — J3 plane from the numerical diagonalization 
data for finite systems by use of the Lanczos algorithm. 



PACS numbers: 75.10.Jm, 75.40. Cx, 75.50.Ee, 75.50.Gg 
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1. Introduction 

In recent years low- dimensional quantum spin systems have attracted a great deal of 
attention. Very recently, Ishii et al. have experimentally studied a trimerized quantum 
spin chain Cu3Cl6(H20)2-2H8C4S02. They have measured the temperature dependence 
of the spin susceptibility, and the magnetization curve at low temperatures. Their results 
show that the magnetic susceptibihty x behaves as x ^ at T ^ and there exists 
the critical magnetic field where the magnetization rises up from zero. Thus they have 
concluded that the ground state of this substance is nonmagnetic. From the result 
of structure-analysis experiment 0, the chain is composed of stacked S = 1/2 Cv?^ 
trimers and separated from other chains by large molecules H8C4SO2. Therefore this 
substance is thought to be well modeled by independent chains of stacked trimers and 
they have proposed a model shown in figure 1. 




Figure 1. Sketch of the model of Cu3Cl6(H20)2 • 2H8C4SO2. Sohd hnes denote the 
intra-trimer couphng Ji, wavy hnes the inter-trimer couphng J2, and dotted hnes the 
inter-trimer couphng J3. 

In this paper we call this model the "distorted diamond (DD) chain model". The 
Hamiltonian of this model is written as 

j j 

+ J3 X! i^3j-2 ■ S^j + S^j ■ S'3j+2) (1) 

j 

where three spins Ssj-i, S^j and /Ssj+i form a trimer. All the coupling constants are 
supposed to be positive (antiferromagnetic). Although it is thought that Ji > J27 -^3 
in Cu3Cl6(H20)2 • 2H8C4SO2 because of its structure, we do not restrict ourselves to 
this case. We note that the point {J2, J3) is equivalent to the point {J2/ J3A/ J3) by 
interchanging the role of Ji and J3. Hereafter we take Ji as the energy unit and set 
J2 = J2/ Ji and J3 = J3/ Ji. 

If we transform the Hamiltonian (P into the fermion representation through 
the Jordan- Wigner transformation, we can see that the fermionic band gap exists at 
M = Ms/3 but not at M = 0, where Ms is the saturation magnetization 0, Thus 
the trimerization itself cannot be the direct reason for the nonmagnetic ground state. 
This can also be explained by considering the necessary condition for the appearance of 
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magnetization plateau proposed by Oshikawa, Yamanaka and Affleck 0, 

n{S — (m)) = integer (2) 

where n is the periodicity of the ground-state wave function, S the magnitude of spins 
and (m) the average magnetization per one spin in the plateau. The periodicity of 
the Hamiltonian (p itself is 3. We see that n = 3 does not satisfy the condition 
(0) with 5 = 1/2 and (m) = 0. Then, if the present model is applicable to 
Cu3Cl6(H20)2 ■ 2H8C4SO2, its ground-state wave function should have the periodicity at 
least n = 6 due to the spontaneous symmetry breaking. In this paper, we explain why 
the non-magnetic ground state is realized and draw the ground-state phase diagram on 
the J2 — J3 plane. 

The Ji = J3 case of the present model was named the "diamond chain" and 
investigated by Takano, Kubo and Sakamoto (TKS) [§. They concluded that the 
ground state of the diamond chain is composed of three phases; the ferrimagnetic phase 
(M = Ms/3) for J2 < 0.909, the tetramer-dimer phase for 0.909 < J2 < 2 and the 
dimer- monomer phase for J2 > 2. The relation between the present model and TKS's 
model will be discussed later. 

This paper is organized as follows. In §2, we explain the mechanism for the non- 
magnetic ground state by use of the analytical method and a physical consideration. 
In §3, we obtain the phase diagram from the numerical data of diagonalization of the 
Hamiltonian (|I]) for finite systems by use of the Lanczos algorithm. The last section is 
devoted to discussions. 

2. Analytical and physical approach 

We consider three special cases at first. When J2 = 1 and J3 = 0, the present model 
is reduced to the simple S = 1/2 chain with nearest-neighbor interactions, the ground 
state of which is the spin-fluid (SF) state, as is well known. In case of J2 = 0, the ground 
state may be ferrimagnetic (M = Ms/3), because the state with S^j =| and ^'3^-1-1 =t 
is the classical ground state. At the point J2 = J3 = 0, the chain is truncated into an 
array of independent trimers. 




Figure 2. Present model in the J2 — 1 case in the linear chain form. 

Next, let us consider the J2 = 1 case. If we re-draw the model in the single chain 
form as in flgure 2, we see this is closely related to the next-nearest-neighbor (NNN) 
interaction model in flgure 3. In fact, the model of flgure 2 is obtained from that of flgure 
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Figure 3. Next-nearest-neighbor (NNN) interaction model. 



3 by removing one NNN interaction of every three NNN interactions. The important 
point is that every spin feels the frustration. 

The NNN interaction model is one of the most important models having the 
frustration and is extensively studied H, ^, |10], 0, [T^, |l^, |l^, |15[. When Jnnn = 
-'nnn/'^i = 0.5, the ground state of the NNN interaction model is an array of independent 
singlet dimers, where the translational symmetry by one spin spacing is spontaneously 
broken ^. Then there should exist the critical point of the ground-state phase 
transition between the spin-fluid (SF) state and the dimer state. Okamoto and Nomura 
numerically determined this SF-dimer critical point, Jnnn = 0.2411. 
Since the model of figure 2 is very similar to that of figure 3, as stated, the SF-dimer 
transition may occur also in the model of figure 2 when J2 is increased. This can be 
confirmed by the bosonization technique in the following way. The effective Hamilton 
of the model of figure 3 in the continuum limit is written as |P, |lT], |12], 0, |T5[ 

2' 



H = — 
2n 



dx 



2tt 



dx cos V2(j) 



(3) 



where Vs is the spin wave velocity and K the quantum parameter which governs the 
algebraic decay of the spin correlation functions 

{S^S^) ~ r-^, {S^S;) ~ r-'/^ (4) 

in the SF state. Due to the isotropic nature of our model, the renormalized value of K 
should be i^' = 1. The variables (f){x) and n(a;) are mutually conjugate, 

[(f){x),U{x')]=i6{x-x') (5) 

The coefficient of the cos-term, in equation (3), is 

y^oc A- 3 Jnnn (6) 

where A is the XX Z anisotropy defined by J^/ J^ which is equal to unity in our isotropic 
model. For the model of figure 2, we can obtain the effective Hamiltonian of the same 
form as (0), but with 

y^cxA- 2 J3. (7) 

We note that the expressions (|^) and (|^ are valid only in the lowest order of A, Jnnn and 
J3. Since we take the continuum limit in the course of deriving the effective Hamiltonian, 
the difference between two models appears as the difference in the expression of y^. We 
note that the spin wave velocity Vs is slightly different between two models, but it does 
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not bring about any essential effect. Thus we can conclude that the model of figure 2 also 
shows the SF-dimer phase transition. Since 3Jnnn in equation (||) is replaced by 2J2 in 
equation (|W), the critical value J2"'' is naively obtained by letting SJ^^n = 2^2'^'^^ which 
leads to J2 — 0.36 by using J^nn = 0.2411. In fact, as shown in §3, the numerical 
result is J2 — 0.354. The dimer ground-state wave function in the model of figure 2 
is two-fold degenerate with periodicity n = 6 and is shown in figure 4. 




Figure 4. Dimer configurations in tlie ground state of the DD chain modei with 
J2 = 1 in the singie chain form. Two spins in an eUipse form a singlet dimer pair 

When J2 7^ 1, we have to take the trimerization effects into the effective 
Hamiltonian. However, in the M = subspace, the trimerization does not bring about 
the mass-generating term such as cos-term in equation (^), although it slightly modifies 
the spin wave velocity Vs- Then, as far as the trimerization is not so large (not so far 
from the J2 = 1 line), the DD chain model also exhibits the SF-dimer phase transition. 
The dimer configuration of the DD chain model is easily known by tracing back of the 
model mapping, which is shown in figure 5. This ground-state wave function is also 
two-fold degenerate with periodicity n = Q. 




Figure 5. Dimer configurations in the ground state of the DD chain model. Two 
spins in an ellipse form a singlet dimer pair. 

Here we summarize the critical properties of the SF-dimer transition using the 
effective Hamiltonian having the sine-Gordon form. The renormalization group 
calculation leads to 
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where L is an infrared cutoff, and 
VqCcK -I 



(9) 



The flow diagram of is shown in figure 6, from which we see that the SF-dimer transition 
is of the Berezinnskii-Kostehtz-Thouless type, as is well known. Since our model is 
isotropic, the renormalized value of K should be equal to unity, as already stated. 
Then, when the system is in the SF state, the starting point of the renormalization 
lies on the SF-Neel boundary line which flows into the origin where K = 1, and moves 
as A — s> B — s> C, as J2 increases. Finally the SF-dimer transition takes place when the 
starting point arrives at the origin. In the SF state, y(j,{L) = —i/q^L) in equation (H), 
resulting in 



yo{L) 



(0) 



,{o)- 



(10) 



y6-'ln(L/Lo) + l 

where i/q'^ is the bare value of yo{L) and Lq is the cutoff length. Then, there appear 
logarithmic corrections in various physical quantities at every place in the SF region 
in our isotropic (i.e., SU{2) symmetric) model. On the SF-dimer critical point (origin 
O), on the other hand, the logarithmic corrections vanish because = y^, = 0. This 
is very peculiar to the isotropic case. Since the SF-dimer transition occurs at y^ = 0, 
one may think that the SF-dimer critical point can be obtained from equations (^ or 
(0). However, as stated, the expressions (H) and (|^) are valid only in the lowest order 
of A, Jnnn and J3, although the critical properties are well expressed by the effective 
Hamiltonian. Then it needs the numerical calculation for determining the SF-dimer 
critical point even in case of J2 = 1. 




Figure 6. The renormalization flow of the effective Hamihonian (^). The phase 
boundaries are shown by thick hnes. In our isotropic case, the starting point of the 
renormahzation hes on the SF-Neel boundary line. As J2 increases, the starting point 
moves as A ^ B ^ C ^ O. If the system has the XXZ symmetry with A = J^' / < 
1, the starting point of the renormalization moves as A' ^ B' ^ C ^ D'. 
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If tfie model has the XX Z symmetry (no longer isotropic) with A = j < 1, 
the starting point of the renormalization moves as A' ^ B' ^ C ^ D' as J2 increases. 
When the starting point arrived at D', the SF-dimer transition occurs. Then, in the 
XXZ symmetric case, the logarithmic corrections exist only at the SF-dimer critical 
point, and do not exist in the SF region. 



3. Numerical result 

To confirm the consideration in §2 and to obtain the ground-state phase diagram on 
the J2 — J3 plane, we performed the numerical diagonalization for finite systems for 
= 6, 12, 18, 24 by use of Lanczos algorithm under the periodic boundary conditions. 
It is very easy to distinguish whether the ground state is ferrimagnetic (M = Ms/3) or 
M = from the numerical data. However, it is difficult to detect the SF-dimer critical 
point from the numerical data, because this transition is of the Berezinskii-Kosterlitz- 
Thouless type |P, with pathological critical behavior. Okamoto and Nomura (ON) 



1^ developed a method by use of which the SF-dimer critical point of the S = 1/2 
NNN interaction model of figure 3 can be successfully determined from the numerical 
data for the energy gaps. Let us explain this method , focusing on its physical meaning. 
In usual cases, the ground state is unique (not twofold degenerate) in finite systems, 
except for the special cases such as the Ising model and the Majumdar-Ghosh model 
1^, 1^. How the twofold degenerate ground state is realized in infinite systems? The 
energy gap of a low-lying excited state of finite systems rapidly decreases as the system 
size increases, and finally degenerate to the ground state in A^ — > 00. Then the linear 
combination of the ground state and the above-mentioned excited state results in the 
twofold degenerate ground state of the infinite systems. In our case, the ground state 
of finite systems has the property S'tot = as far as it lies in the M = subspace 
(i.e., except for the ferrimagnetic case). The twofold degenerate dimer state of infinite 
systems also has the property S'tot = 0. Then the above-mentioned excited state should 
be also of S'tot = 0, because of the law of the addition of the angular momentum. Then 
we can conclude that the lowest excitation in finite systems is of S'tot = in the dimer 
region. In the SF region, on the other hand, the lowest excitation should be of the 
spin-wave type with S^^^ = ±1 (one magnon state). In the present case, the excitation 
with the same energy exists in the St^t = subspace due to the isotropic nature. This 
means that the system has three- fold degenerate lowest excitation with S'tot = 1, when 
it lies in the SF region. 

From the above physical consideration, we can write down the criterion 

AEss(A^) < AEst(A^) dimer state 

AE,,{N) > AEst(A^) spin fiuid state (11) 

where Aii^ss(N) and Ai?st(N) are the singlet-singlet energy gap and singlet-triplet energy 
gap for a finite-size system with A^ spins, defined by 

AE,,{N) = Ei{N, S(*°*) = 0) - Eg{N) (12) 
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A£;,,(iV) = Eo(iV, S^'°'^ = 1) - Eg(iV) (13) 

respectively. Here Eq{N,S^^°^^) and Ei(A^, are the lowest and second lowest 

energies in the subspace with 5*'-*°'^'' and Eg = £'o(S'*^*°*^ = 0), respectively. This criterion 
can be obtained also by use of the effective Hamiltonian representation, renormalization 



group method and the conformal field theory |12, 14, 15 1. 



I 
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Figure 7. Level crossing between AEss{N) and AEss{N) when J2 = 0.8 and N = 18. 
From the crossing point, we obtain J^"\n — 18) — 0.375263. 




~(C1') ~ 

Figure 8. The extrapolation of J3 to TV ^ cx) in J2 = 0.8 case. From this, we see 
J^"^ = 0.352 ±0.001. 

Figure 7 shows the crossing between S'tot = and S'tot = 1 excitations when 
J2 = 0.8 and = 18. By use of the interpolation, we see that the crossing point 
is J^^'^\n = 18) = 0.375263. We can obtain the SF-dimer critical point of the infinite 
system by extrapolating J^^^\n) to — > 00. The finite-size dependence of J^^^\n) has 
the form 

jt\N) = jt\oo) + (const /iV2) (14) 
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8 shows the extrapolation of j'^^\n) to 



due to the existence of the irrelevant fields, as was discussed in 0, |T^. Figure 

oo in case of J2 = 0.8, resulting in 



Jl' = 0.352 ± 0.001. 



By sweeping parameters, we finally obtain the phase diagram on the J2 — J3 plane. 
The result is shown in figure 9. We note that the point {J2, J3) is equivalent to the point 
{J2/J3, ^/ J3) by interchanging the role of Ji and J3, as already stated in §1. 




Figure 9. Phase diagram of the DD chain model. The estimated errors in the critical 
values are less than 0.001. The J3 = 1 case is reduced to the model of Takano et al. 
(see §4). 



4. Discussion 

In §2, we have stated that the ground-state quantum phase transition of the DD chain 
model (the present model) has the same universality class as that of the NNN model 
(figure 3). We have confirmed this analytically by use of the effective Hamiltonian 
representation. Here we also confirm this by numerical method. 




Figure 10. The interpolation model between the DD chain model with J2 = 1 and 
the NNN interaction model. Solid lines denote Ji, dotted lines J3 and dot-dashed lines 
J4. 

Let us consider the model of figure 10 which interpolates between the DD chain 
model with J2 = 1 (figure 2) and the NNN interaction model (figure 3). When J4/ J3 = 
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the interpolation model is reduced to the DD chain model with J2 = 1, and when 
Ji/ Jz = 1 to the NNN interaction model. Figure 11 shows the level crossings between 
Ast(iV) and Ass(A^) for J4/ J3 = and J4/ J3 = 1 cases. The behavior of level crossing is 
essentially the same when the parameter J4/ J3 runs from J4/ J3 = to J4/ J3 = 1, as can 
be seen from figure 11. Furthermore any other excitations cross them between J4/ J3 = 
to J4/J3 = 1. Figure 12 shows the SF-dimer critical point J^^ of the interpolation 
model, in which the critical point smoothly changes. Thus we can safely conclude that 
the ground-state quantum phase transition of the DD chain model (the present model) 
has the same universality class as that of the NNN interaction model (figure 3). The 
DD chain model with J2 = 1 is obtained from the NNN interaction model by removing 
one NNN interaction in every three NNN interactions, as already stated. Instead of 
removing, a similar (but not exact) effect may be realized by decreasing the strength of 
the NNN interaction to 2/3 of the original strength. If this is the case, the SF-dimer 
critical point of the DD chain model with J2 = 1 is 3/2 of that of the NNN interaction 
model, which results in J^^^ = (3/2) x 0.2411 ~ 0.36. This semiqualitatively explains 



our numerical result J^^^^ 
equations (H) and (^. 



0.354 ± 0.001 when J2 = 1. This fact also appears in 



■ooo 
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oo 
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2 
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0.4 



Figure 11. Level crossing of the interpolation model when J4/J3 = and J4/J3 = 1 
for N = 18. Closed circles and squares represent Ast(iV), and open circles and squares 
As.(iV). 

Takano, Kubo and Sakamoto (TKS) [H investigated the J3 = Ji case of the present 
DD chain model (see figure 13(a)). They concluded that the ground state of their 
model is composed of three phases. The ferrimagnetic phase (M = Ms/3) appears when 
J2 < 0.909. In the tetramer-dimer (TD) phase, which appears when 0.909 < J2 < 2, 
the state is exactly the regular array of tetramers and dimers as shown in figure 13(b). 
Figure 13(c) shows the dimer-monomer (DM) state, appearing when J2 > 2, which is 
composed of the regular array of the singlet dimers and free spins. Because of the free 
spins, the DM state is macroscopically degenerate. 
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Figure 12. The SF-dimer critical points of the interpolation model between the DD 
chain model with J2 = 1 and the NNN interaction model. 



Let US discuss the relation between our model and TKS's model. In the DM state of 
TKS's model, the monomers are completely free but this is very peculiar to this model. 
In our model, since the symmetry of a diamond is broken because J3 7^ 1, the monomer 
is no longer free and has an effective interaction between neighboring monomers through 
the dimer between them. Therefore the DM state of TKS's model is smoothly connected 
to the spin- fluid state of our model, as can be seen in figure 9. The tetramer in the TD 
state is also special to TKS's model. When the symmetry of the diamond is broken, the 
tetramer is decomposed into two dimers existing on stronger bonds, as is shown in figure 
5. Then the TD state of TKS's model is a special case of the dimer state of our DD 
chain model. Thus the physical pictures of our model and TKS's model are consistent 
with each other. 

We can confirm the reliability of our numerical results by checking the properties 
of excitations. At the SF-dimer critical point, the cos-term of the effective Hamiltonian 
(I) vanishes and the relation A£'ss(A^) = AEst{N) holds, as discussed by Okamoto 
and Nomura 0, 0, |15|. Then the system is purely Gaussian and has the low-lying 
excitation energies proportional to in finite systems. The lowest order correction 
to 1/A^ may be of This correction comes from the band curvature and the wave- 

number dependence of the coupling constant of the interaction between Jordan- Wigner 
fermions, which was neglected in the course of deriving the effective Hamiltonian. Figure 
14 shows the l/A^^-dependence of NAE{N) at J2 = 1.2 and J3 = 0.39, which is the 
SF-dimer critical point where Aii^ss(oo) = Ai?st(oo). As can be seen from figure 14, the 
size dependence of the lowest excitation is well expressed as 

NAE{N) =a + (b/N^) (15) 

which is consistent with the above-mentioned discussion. The quantity a is related to 
the spin wave velocity as 



(16) 
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(b) 




c 




Figure 13. (a) The model of Takano, Kubo and Sakamoto, (b) The tetramer-dimer 
(TD) state. The rectangles represent tetramers and the ellipses singlet dimers. (c) 
The dimer-monomer (DM) state. 



where x is the scahng dimension of this excitation, which is equal to 1/2 at the critical 
point WA III, Since a = 3.851 ± 0.001 in case of figure 14, we obtain 



1.226 ±0.001 



;i7) 



3.90 



3.85 




0.01 



Figure 14. System-size dependence of the scaled excitation gap NAE on the SF- 
dimer critical point. 



The system-size dependence of the ground-state energy also provides us with useful 
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information. Under periodic boundary conditions, it is written as |T6| , |T7 

EJN) , , TTfsC , , 

where Eg{N) is the ground-state energy of the A^-spin systems, eg(oo) the ground state 
energy of the infinite system per spin, Vs the spin wave velocity, and c the conformal 
charge which is equal to unity in our universality class. Figure 15 shows the system-size 
dependence of the ground state energy in case of J2 = 1.2 and J3 = 0.39. From the 
slope of the line, we obtain 

fs = 1.23 ±0.01 (19) 

which well agrees with equation (|T^. The fact that the values of Vs given by equations 
(p!7|) and (|T9|) agree within the numerical error justifies our numerical analysis. 



-0.425 



5 



-0.430 



-0.435 




1/N 



0.01 



Figure 15. System-size dependence of the ground-state energy Eg{N) on tiie SF- 
dimer critical point. 

The spin wave velocity Vs can be also obtained from the lowest excitation having 
S^^^ = and A; = 2ti/N by 

_ J. ^ iVAE(iV.54. = 0,t-^2ViV) po) 

We note that this formula is free from the logarithmic corrections even in the SF 
region [|T^. Figure 16 shows this extrapolation procedure, which brings about Vs = 
1.228 ± 0.001. This is also consistent with equation ( |17|) and equation (pJl) . 

Let us discuss the logarithmic corrections in the SF region. As stated in §2, there 
appear logarithmic corrections in various physical quantities when the system is in the 
SF region. This is very peculiar to our isotropic case. In the following we check this 
point numerically. The singlet- singlet gap and the singlet-triplet gap are expressed as 

AE^N) = AE^m = (21) 

respectively, where Xss and the scaling dimensions 

^ss = ^ (1 + IvoiN)^ ^st = ^ (1 - IvoiN)'^ (22) 
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Figure 16. Extrapolation procedure of equation (pOj). From the intersection, we 
obtain = 1.228 ± 0.001. 



with yo{N) given in equation (0). It is difficult to directly detect the logarithmic 
dependence in equations (^TJ) and (|22|) from the numerical data for AEss{N) and 
AEst{N). Because the logarithmic corrections are very slowly varying with respect 
to the system size A^, its effects are actually observed as the change in the spin wave 
velocity Vs between AEgg and AEst- As an example, let us take the J2 = 1.2 and 
J3 = 0.35 point which lies in the SF region. In fact, as shown in figure 17, the spin wave 
velocities are estimated to be Vs = 1.362 ± 0.001 and Vs = 1.248 ± 0.001 from AEss and 



AEst, respectively. Okamoto and Nomura |12, 04, used the "averaged excitation 



AE,,,{N) = - {AEss{N) + 3AEgt{N)} 



(23) 



in which the lowest order logarithmic corrections vanish, as can be seen from equation 
(p^). From figure 17 we obtain Vg = 1.277 ± 0.001 by use of A_Eavc- When we calculate 
the spin wave velocity through the formula (^), we obtain Vg = 1.277 ± 0.001, which 
shows very good agreement with that from Ai^ave- We note that there is no logarithmic 
correction in formula (^), as already stated. Thus our numerical analysis is consistent 
with our consideration in §2 with respect to the logarithmic corrections, although we 
could not directly observe its existence. It may need systems with several thousand (or 
more) spins to directly observe the contribution of the logarithmic corrections. 

In summary, we have explained that the frustration brings about the non-magnetic 
ground state in our DD chain model by use of the analytical method, physical 
consideration and numerical method. We have also obtained the phase diagram on 
the J2 — J3 plane numerically. 
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Figure 17. Behavior of AE'ss, AEgt and Ai?ave at the point (J2, ./a) = (1-2, 0.35) in the 
SF region. From the intersection of the ASave hne, we obtain Vg = 1.217 ± 0.001. The 
apparent spin wave velocities are = 1.362 ± 0.001 from AEgg and Vg = 1.248 ± 0.001 
from AiJgf 



the numerical calculation was done by use of program package TITPACK Ver.2 coded 
by H. Nishimori. 
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